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Effects of Downscattering on the Continuum and Line Spectra in 
Powerful Wind Environment. Monte Carlo Simulations, 
Analytical Results and Data Analysis 

Philippe Laurent^'^ and Lev Titarchuk'^'^ 
ABSTRACT 

In Paper by Titarchuk & Shrader the general formulation and results for 
photon reprocessing (downscattering) that included recoil and Comptonization 
effects due to divergence of the flow were presented. Here we show the Monte 
Carlo (MC) simulated continuum and line spectra. We also provide an analytical 
description of the simulated continuum spectra using the diffusion approximation. 
We have simulated the propagation of monochromatic and continuum photons 
in a bulk outflow from a compact object. Electron scattering of the photons 
within the expanding flow leads to a decrease of their energy which is of first 
order in V/c (where V is the outflow velocity). The downscattering effect of 
first order in V/c in the diverging flow is explained by semi-analytical calcula- 
tions and confirmed by MC simulations. We conclude that redshifted lines and 
downscattering bumps are intrinsic properties of the powerful outflows for which 
Thomson optical depth is greater than one. We fitted our model line profiles to 
the observations using four free parameters, P = V/c, optical depth of the wind r, 
the wind temperature kT^, and the original line photon energy Eq. We show how 
the primary spectrum emitted close to the black hole is modified by reprocessing 
in the warm wind. In the framework of the our wind model the fluorescent iron 
line Kq is formed in the partly ionized wind as a result of illumination by central 
source continuum photons. The demonstrated application of our outflow model 
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to the XMM observations of MCG 6-30-15, and to the ASCA observations of 
GRO J1655-40, points out a potential powerful spectral diagnostic for probes of 
the outflow-central object connection in Galactic and extragalactic BH sources. 

Subject headings: accretion — stars: radiation mechanisms: nonthermal 

1. Introduction 

Titarchuk & Shrader (2005), Paper I, hereafter TSh05, have shown that downscattering 
modification of the primary photon spectrum by an outflowing plasma is a possible mecha- 
nism for producing the continuum excess in the ~ 10 keV spectral region. This continuum 
excess can be formed only if the energy index of the primary photon spectrum is less than one 
(see also Sunyaev & Titarchuk 1980, hereafter ST80). This is usually attributed to Comp- 
tonization by a static reflector, such as a downward- or obliquely-illuminated accretion disk, 
although the overall continuum form differs from that of the basic Compton reflection form. 
TSh05 thus suggest, that in at least some cases, the outflow downscattering effect rather 
than the standard Compton reflection mechanism is responsible for the observed "excess" 
hard- X-ray continuum. 

To obtain a description of the X-ray photon spectrum, many authors (see e.g. Gilfanov, 
Churazov Sz Revnivtsev 1999, hereafter GCR99, Pottschmidt et al. 2003, and etc) used an 
empirical model in which each source spectrum is a sum of power law spectrum with photon 
index F and a multi-temperature disk blackbody. To this continuum, a reflection spectrum 
after Magdziarz & Zdziarski (1995) was added. However, Shaposhnikov & Titarchuk (2006), 
hereafter ST06 argue that these spectral features arc results of reprocessing of the central 
source hard radiation in the outflow (warm wind) that surrounds the central BH. This claim 
is supported by observed power decay of quasiperiodic oscillations (QPOs) towards the softer 
states that is presumably a result of reprocessing of time signal in the extended relatively 
warm wind (see Fig. 6 in ST06). 

The basic idea is that electron scattering of photons from a central source entering the 
expanding outflow experience a decrease in energy (downscattering) . The magnitude of this 
decrease is of flrst order in V/c and in E/nieC^ where V is the outflow speed, c is the speed 
of light, E is the initial photon energy and mg is the electron rest mass. 

TSh05 have developed an analytic formulation for the emergent spectrum resulting 
from photon diffusion in a spherically expanding Comptonizing media, characterized by two 
model parameters: an average number of scatterings in the medium Nav, and efficiency the 
energy loss in the divergent flow e. In this formulation of the Radiative Transfer problem. 
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the number of photons emitted in the wind shell equals that which escape to the Earth 
observer (i.e. the photon number is conserved). Consequently, if the high energy photons 
lose their energy on the way out but the number of photons is conserved, it has to result in 
the accumulation of the photons at a particular lower energy band. 

The basic idea of the wind photon downscattering is depicted in Figure 1. There presents 
a simple explanation of the diverging flow effect on the photon propagation through the 
medium. A photon emitted outwards near inner boundary and then scattered at a certain 
point by an electron moving with velocity Vi, is received by an electron moving with velocity 
V2 as shown with frequency V2 = '^i [1 + (Vi — V2) ■ n/c] where n is a unit vector along the 
path of the photon at the scattering point. In a diverging flow (Vi — V2) ■ n/c < and 
photons are successively redshifted, until scattered to an observer at infinity. The color of 
photon path (in Figure 1) indicates the frequency shift in the rest frame of the receiver 
(electron or the Earth observer). On the other hand, referring to the right-hand side of 
Figure 1, in a converging fiow (Vi — V2) • n/c > and photons are blueshifted. 

In this Paper we present the results of extensive Monte Carlo simulations of the down- 
scattering effects in the wind (diverging fiow). In addition to the scattering, Doppler and 
Compton effects we include a photoabsorption effect in the diverging fiow (wind) . Any pho- 
ton with energy of about 7-8 kcV and higher interacting with outflow plasma is more likely 
to be absorbed by the flow and be reemittcd at energies about 6.4-6.6 keV, depending on the 
ionization stage of the flow (sec Kallman ct al. 2004). Laming and Titarchuk (2004), here- 
after LT04, calculated the temperature and ionization balance in the outflow determined by 
seeking the temperature at which the outflowing gas attains photoionization-recombination 
equilibrium. The gas is heated by Compton scattering and photoionizations by photons from 
the central compact object and is cooled by radiation, ionization, and adiabatic expansion 
losses. The solution of this problem is general regarding the calculation of the temperature 
and ionization balance in the target illuminated by the hard radiation. It potentially has 
many applications beyond Astrophysics (for example, in Radiation Safety Physics). Our 
simulations allow us to reveal the spectral, timing and absorption properties of the emergent 
radiation. 

We update the analytical results of the downscattering modifications of the continuum 
obtained by TSh05 (using Fokkcr-Planck technique) and we verify them using the Monte 
Carlo (MC) simulated spectra. We demonstrate that the analytical formulas for the contin- 
uum spectra well describe the MC spectra almost for any j3 = V/ctO.8. 

Many authors [e.g. Tanaka et al. (1995); Nandra et al. (1997); Wilms et al. (2001); 
Uttley et al. (2004); Miller ct al. (2004a,b)] have found unusual curvature (red-skewed 
features) at energies between 2 and 8 keV in X-ray spectra of a number of galactic and 
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extragalactic BH sources. Significantly the fluxes of these features remain nearly constant 
despite the large changes in the continuum flux as shown in the observations of MCG-6-30- 
15, NGC 4051 and others [Markovitz, Edelson, & Vaughan (2003)]. This lack of variability 
means that fast changes (which are almost certain to occur near the BH) are not obtained. 
These observations motivate us to make a detailed study of the spectral line formation in 
the outflows that possibly occur in black hole sources. In this Paper we demonstrate that 
red-skewed features can be a result of scattering of iron line photons in the powerful outflow. 

We offer a model in which the hard radiation of the central object illuminates and heats 
the outflow (wind) region originated in the outskirts of the disk (well outside the innermost 
part of the accretion disk near the BH). The line photons are generated and scattered in the 
outflow. Our basic scenario is illustrated in Figure 1. The wind originates at a distance r^n 
from the central black hole and is of a density such as to give a Thomson scattering optical 
depth tt close to unity far from the black hole. The optical depth in the Fe K continuum 
is about 1-3 times higher than that due to electron scattering (assuming a solar abundance 
of Fe, and depending on charge state, see Kallman et al. 2004), and so Fe Ka formed by 
inner shell ionization of Fe ions in the outflow by the continuum from the central black hole 
only comes from a smaller inner region. In fact, using the calculations of the ionization state 
of the outflow shell by LT04, we find that K^— photon sources are exponentially distributed 
over a skin at the bottom of the outflow, the skin thickness being 4% of the optical depth 
of the wind shell. Furthermore because the iron photo-absorption opacity ax is about few 
times the Thomson opacity, (Tt, any photon around the K-absorption edge energy must be 
absorbed. In fact, these edges are seen with the strong Kq, hues during X-ray superbursts 
detected from NS source 4U 1820-30 (see Figs. 5 and 9 in Strohmayer & Brown 2002). 

A thorough analysis and review of the diffusion theory of photon propagation in an 
optically thick fluid in bulk outflow has been provided by Titarchuk, Kazanas & Becker 
(2003), hereafter TKB. They show that the iron hue is produced in an effectively optically 
thick medium. Its red wing is the result of multiple scattering, each scattering producing a 
first order y/c— redshift. This process produces a red wing to the line without any particular 
fine-tuned geometric arrangement. 

The TSh05 and TBK results are obtained using the Fokker-Planck equation [Blandford 
& Payne (1981) see also Lifshitz & Pitaevsky (1981) for a general treatment of the particle 
acceleration problem] which is valid for (3 = V/c 1 and Tg// > 1. The MC method does 
not have these limitations. 

We show the results of MC simulations of scatterings of continuum and monochromatic 
photons for various values of plasma temperatures kTe, outflow bulk velocity f3, and Thomson 
optical depth r. The details of the MC method applied to the outflow model are given in §2. 
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We compare the results of the MC simulated continuum with their analytical description and 
discuss the downscattering modifications of the incident spectrum in §§3.1-3.2 The results of 
the simulated line formation in the outflow and its explanation using the analytical theory 
are shown in §3.3. Examples of the model fits to the XMM data obtained for MCG-6-30-15 
and to ASCA archive data of galactic BHs are shown in §4. The results of the simulations 
of overall spectra including the scattering, photon absorption effects are presented in §5. 
Discussion and conclusions follow in §6. 

2. Details of the Monte Carlo Method 

Laurent & Titarchuk (1999), hereafter LT99, simulated the Comptonization of photons 
using a Monte-Carlo method in the Schwarzschild metric. The geometry used in these 
simulations consists of a spherical shell harboring a black hole in its center. Once given the 
cloud Thomson optical depth t, is an electron number density measured in the local rest 
frame of the flow that is 

= r{rs/rf/[rs(rT{rs/rin - rs/w)]- (1) 

Here ar is the Thomson cross section. Tin and Tout are the inner and outer radii of the shell 
respectively, and rs is the Schwarzschild radius. Simulations made for different values of rj„ 
and Tout (that are much greater than the Schwarzschild radius) have shown no change of 
the line proflle and continuum provided the shell optical depth remains the same. We have 
varied the shell electron temperature from 0.01 to 1 keV. Along with the thermal motion of 
the electrons in the cloud, we have also taken into account the bulk velocity of the outflow 
(5 — V/c. In these simulations, the spherical wind has a constant bulk velocity within the 
flow. It is worth noting that an increase of the shell temperature from 0.01 to 1 keV slightly 
increases the blue wing of the observed line, but induces no change in the red wing that we 
are mostly interested in here. 

For the timing studies described in the next section we have also computed the path 
length that each photon in the simulations has covered until it escapes from the shell. 

3. Dowscattering effect in an outflow 

The problem of photon propagation in a fluid in bulk motion has been studied in detail 
in a number of papers [see e.g. Blandford & Payne (1981), Payne & Blandford (1981), 
Nobih, TuroUa & Zampieri (1993); Titarchuk, Mastichiadis & Kylafls (1997); TKB; TSh05; 
Laurent & Titarchuk (1999), (2001). In particular, TSh05 present a general formulation and 
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a solution of the spectral formation in the diverging outflow. They demonstrated that the 
resulting spectrum can be formed as a convolution of energy and spatial diffusion solutions. 

To apply these Radiative Transfer results to observational data, TSh05 have developed a 
generic analytical formulation which leads to a simple analytic expression for modification of 
the emergent spectrum due to recoil and velocity divergence effects in the flow. For the recoil 
effect, TSh05 extend the results by Sunyaev & Titarchuk (1980), hereafter ST80, who show 
that the downscattering feature, or "bump", can appear superposed on power-law spectra 
with energy spectral indices a < 1 as a result of photon diffusion through a static cloud. 
From ST80 and calculations presented in TSh05, one might conclude such bumps are not a 
necessarily a feature of disk reflection, but they can be also formed as a result of the photon 
reprocessing in a relatively cool, ambient plasma characterized by temperatures of order 10^ 
K. 



3.1. Time distributions of escape photons 

In order to understand and explain the shapes and the main features of the MC simulated 
spectra we attempted to reproduce them using the analytical diffusion approximation. In 
framework of the diffusion theory (TSh05) the resulting spectrum can be presented as a 
particular convolution of a primary photon spectrum ip{z) and a time distribution of escape 
photons V{u) (see Eq. |2l below). Here u = t/tfp, t and tfp = l/c are the time for photon to 
escape from the wind and the mean free path time between two consequent photon scatterings 
respectively; z = E/me(? is the dimensionless photon energy. ST80 and Sunyaev & Titarchuk 
(1985) hereafter ST85, demonstrate that the photon escape distribution V{u) exponentially 
decays with m as m > 1 for any bounded configuration, namely V{u) ~ 6exp(— 6m). 

In Figure 2 we present the simulated time distribution V{t/tfp) for optical depth of 
the wind r equal to 4. As it was found by ST80 and ST85, the time distribution decay 
exponentially exp{—bt/tfp) for t > tfp when the photons propagate in any bounded medium. 
This result is used for derivation of the analytical spectra (see Eqs|2]). It means that the best- 
fit b— parameters obtained by fitting exp{—bt/tfp) to the time distributions and those b— 
parameters obtained by fitting the analytical spectra to the simulated spectra are identical 
within errors of b. We computed the 6— values obtained by these ways and we found they 
are identical. This consistency check gives us more confidence of accuracy of our simulations 
and their analytical description. 

In Figure 3 we show the photon distribution for emergent photons over a range of 
deflecting angles with respect of the direction of the incident photons. The outflow shell 
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is radially illuminated be the central source radiation. The simulations are presented for 
r = 1, 2, 4 and (3 = 0.1. In fact, we find in our simulations that the distribution of the 
scattering component weakly depends on r. Practically, there is little difference between that 
related to r = 4, and r = 1, 2. There is only the difference in the normalization of the direct 
(non-scattering) component for which the ratio of An{t = 4:)/An{t) = exp(— 4-|-r). The 
distribution has a strong maximum about 50°. Thus we expect that our results for spherical 
wind can be applicable for the conical wind if the half-angle of the conical wind is about 50° 
or more (wide open jet or wind). 

As we demonstrate in Figure 2 the simulated escape time distribution in the wind has 
an exponential profile which is a typical signature of the photon diffusion propagation in 
any bounded configuration (ST85). This leads us to conclusion that we can also treat the 
radiative transfer in the wind in the framework of the diffusion model (see TSh05). 

3.2. Monte Carlo simulated spectra of the continuum and their analytical 

description 

For the completeness of the present work we reproduce some necessary details of the 
analytical diffusion model by TSh05. Particularly, the downscattering modification of the 
incident spectrum (p{E) can be obtained using formula (24) in TSh05: 

1 rUmax.e{z) 

TAz.e) = - / e-''^/^^'-'\z,uM^'^'\z,u)]V{u)du, (2) 

^ Jo 

where V{u) is the distribution of the photons over their dimensionless escape time u = t/tfp, 

) = (3/e)ln(l + £/32) > 1, (3) 

is a maximal number of scatterings related to the dimensionless energy z = E /m^c^ . 

zo = i^^'\z,u) = (£/3)/[(l + £/3z)exp(-£M/3) - 1], (4) 

is the dimensionless photon energy zq = EQ/nieC^ at m = and e <^ 1 is a coefficient of the 
diverging term in the kinetic (Fokker-Planck) equation (see Eqs. 16, 17 in TSh05 and related 
discussion). In fact, e is an average relative energy redshift per scattering. Below (see Eq. 
[12]), we evaluate these values and verify them using results of our Monte Carlo simulations. 
We find that e ~ /3/2r is of order of 10~^ for (3^0.1 and for r of a few. 

This result is a generalization of ST80 and TKB results. ST80 derived the spec- 
tra when the downscattering effects due to the recoil was taken into account. It should 
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be noted that the downscattering modification of the spectrum occurs when photons un- 
dergo multiple scatterings u. Thus we can expand the integrand function W^^\z,u) = 
e~'^^'^^^[tp^'^\z,u)](p[ip^'^\z,u)] over u in formula ([2]): 

W^'\z, u) ^ W'^'\z, 0) + Wi'^'iz, 0)u + Wi'^^'^z, 0)mV2 + Wi''>"'{z, 0)uV3! + ... (5) 



For the simplest case of the power-law incident spectrum and neglecting the diverging 
term in the kinetic (Fokker-Planck) equation {e = 0) the resulting spectrum (Eq. [2]) is as 
follows 

jr^(^) = - [ (i/z - u)''-^Viu)du. (6) 



z 

Using the expansion ([5]) one can obtain (see Appendix A, Eq. IA6P 

TeU) = Az-'[M„(z-') + Il^,M,(.-i) + '•^ ~ ~ "K Hk(z-') + ... 

+ -"' .'AM.-') + ...]. (7) 

where 

Mn{z'^) = b / u"" exp{-bu)du for n = 0, 1, 2... (8) 
Jo 

Note, the first term in the parenthesis of the right hand side Mo(z~^) = 1 for z <^ 1 and 
the second term of that (1 — a)zMi{z~^) is positive if energy index of the primary photon 
spectrum a < 1. Thus one can expect the photon accumulation bump at low energies for the 
downscattered power-law spectra for which the primary photon spectrum is quite hard (the 
energy index a < 1). Indeed, in the pure scattering case the photon number is conserved, 
so if large numbers of the downscattered high energy photons are removed from the high 
energy part of the incident spectrum, then they should are seen at lower energies as a so 
called downscattering bump. 

In Figures 4, 5 we present the Monte Carlo simulated spectra along with the analytical 
spectra (Eq. [7j) for the power-law incident spectra for r = 4 and r = 2 respectively. 

The downscattering (accumulation) bump and softening of spectrum at high energies 
are clearly seen in these spectra. As optical depth increases more prominent bumps are 
formed in the outflow. The softening of the analytic spectra due to the downscattering 
effect is well described by exponential decay of moments of the time distribution, Mi{l/z) 
{i = 0, 1, 2, ..) at higher energies (i.e. when b/z > 1), see Eqs. (IHl), (1A9I1 . The parameter b is 
only a free parameter to fit the analytical formula, we fix e = for which Umax = 1/^, see 
Eq. ([3]), namely limg^o ^ma^' = 1/^- For r = 4 the analytical spectra ([7]) perfectly represent 
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the simulated spectra for a wide range of the wind velocities V = (3c^0.8c. It is worth noting 
that we need not more than 6 terms of asymptotic series ([7]) to achieve the highest accuracy 
of the analytical spectra calculations. 

For/3 = 0.05, 0.1, 0.3, 0.8 the best-fit parameters 6 = 0.07, 0.07, 0.11, 0.45 respectively. 
The inverse of 6— parameter, Nav = 1/6 is an average number of scatterings undergone by 
the photons in the wind. For a faster wind the number of scatterings naturally decreases 
because photons are effectively taken by the outflow and less photons are scattered off outflow 
electrons. ST85, calculate b-parameters (and the average number of scatterings Nav) for 
spherical and plane (disk) geometry (Table 1, in ST85, where /3 corresponds to our value b). 
ST85 present their results for the static case i.e the zero bulk velocity case. Thus we should 
compare their results with our results when (3 = 0.05, 0.1. Moreover, our wind extended 
envelope of optical depth r is equivalent to the disk of half-optical depth r/2. For example, 
for r = 4, our b = 0.07 should be related to 0.101 in ST85 for r = 2. As expected our b-value 
is smaller (or 1/6 greater) than that in ST85, because the empty internal cavity of the wind 
envelope essentially amplifies the number of photon scatterings with respect to that in a 
plane disk of the same optical depth. The scattering effect of the cavity is suppressed when 
the wind velocity (or /3) increases. For f3 = 0.3, the cavity effect and the wind velocity on 
scattering compensate each other. Our value of 6 = 0.11 is very close to that in ST85. 

In Figure 5 we show the simulated spectra along with the analytical spectra for r = 2. 
One can notice that the best-fit analytical spectra slightly deviate from the simulated spectra 
at high energies. For /? = 0.05, 0.1, 0.3, 0.8 the best-fit parameters 6 = 0.13, 0.15, 0.21, 0.62 
respectively. ST85 give the corresponding value of 6 = 0.234 (related to our 6-parameter for 
f3 = 0.05). From this comparison it is evident that the scattering effect of the cavity becomes 
stronger with the optical depth decrease. As expected, the analytical spectra as solutions 
of the diffusion problem in Thomson regime may depart from the Monte Carlo simulated 
spectra for r 2 at high energies where the electron opacity is less than the Thomson opacity. 
The central source spectrum (dashed line) is slightly less modified by the scattering in the 
wind (see MC histogram for the simulated spectrum) than that using the diffusion spectrum 
(solid line). 

Another interesting case is the downscattering modification of Comptonization spectrum 
(see e.g. ST80 and Titarchuk 1994). The incident Comptonization spectrum can be well 
fitted by a power-law spectrum with an exponential cutoff, namely 

Vcomp{z) ~ Az~°'exp{-z/z^). (9) 

where = E^/nieC^ is the dimensionless photon energy and the cutoff energy i?* is related 
to the Compton cloud electron temperature kT^, namely ^ 2kTe. The downscattering 
modification of the spectrum, (see Eq. IA8|) can be written using the expansion ([5]). It can be 
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shown that for e = and ^ oo formula flA8l) is reduced to formula ([7]). In Figures 6 and 
7 we present the Monte Carlo simulated spectra for the power-law with exponential-cutoff 
incident spectra for r = 4, r = 2 respectively along with the analytically calculated spectra 
(see Eqs. IA8[ IA20tiA25"|) . The MC simulated spectra are well described by the analytical 
model for all f3t0.8. It is worth noting that the downscattering bumps are not seen in either 
the simulated nor model spectra. 

Forr = 4and/5 = 0.1, 0.3, 0.5, 0.8 the best-fit parameters are 6 = 0.06, 0.10, 0.18, 0.38 
respectively. For r = 2 and the same (3 the best-fit parameters are b = 0.14, 0.21, 0.30, 0.59 
respectively. It is worth noting that the values of 6— parameter for small (3t0.1 is almost 
independent of the shape of the incident spectrum (within the error bar Sb = ±0.005) and 
it depends on r only. 

Comparison of the simulated and analytically calculated spectra leads us to conclude 
that the simulated spectra can be well described by the analytical model with only one free 
parameter b that is related to the mean number of photon scattering in the wind (A^^^ = 1/b). 
It means that the continuum spectral formation in the pure scattering wind is dictated by 
the mean number of scatterings Nav only but that the value of Nav is determined by the 
combined effect of the wind optical depth r and the dimensionless wind velocity jS. 

3.3. Kq,— Line formation in the outflow 

Another observational feature of the wind is the strong broad feature of Kq, line in the 
spectrum. Recently, Shaposhnikov & Titarchuk (2006) demonstrate that these strong iron 
lines are present in all spectral states of Cyg X-1. Also they show that the equivalent width 
(EW) of the Ka line increases with the photon index F from about 150 eV in the low/hard 
state to about 1.3 keV in the high/soft and very soft states. 

The X-ray photons of the central source illuminates the wind (see Fig. 1). The wind gas 
is heated by Compton scattering and photoionizations by photons coming from the central 
object and is cooled by radiation, ionization, and adiabatic expansion losses (LT04). The 
photons above the K-edge energy are absorbed and ionize iron atoms which leads to the 
formation of the strong K^^ line. LT04 calculated ionization, temperature structure and the 
equivalent widths of Fe Kq, line formed in the wind. For the wide set of parameters of the 
wind (velocity, the Thomson optical depth r) and the incident Comptonization spectrum 
(the index and the Compton cloud electron temperature) they established that EW of the 
line should be about 1 keV and less for the line to be observed. LT04 also predicted that 
in this case the inner radius of the wind should be situated at (10'^ — 10'^)/ro Schwarzschild 
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radii away from the central object. 

3.3.1. The spectral line redshift in outflow. Analytical Description 

In Figure 1 we show the picture of the hne photon propagation in the wind. At each con- 
secutive scattering the hne photon in average loses its energy, in other words it is redshifted. 
A photon emitted near the inner boundary and subsequently scattered by an electron moving 
with velocity Vi, is received by an electron moving with velocity V2 as shown in Fig. 1. 
The change in frequency is 

1^2 = i^i[l + (Vi - V2) • n/c] (10) 

where n is a unit vector along the path of the photon scattered at the next point. 

We remind a reader that in the static and spherical-symmetric medium and where elec- 
trons experience a random (Brownian) motion only, the mean energy change per scattering, 

< AE >, for photons undergoing quite a few scatterings is proportional to square of electron 
velocity, namely < AE >oc (V/c)'^ (see e.g. Rybicki & Lightman 1979 and ST80). 

On the other hand for photons undergoing numerous scatterings in the bulk flow < 
AE > is just proportional to V/c. [see Titarchuk, Mastihiadis & Kylafls 1997, Appendix D 
and TKB). In order to make estimate of the loss per scattering, one has to flnd the response 
of the energy operator of the diffusion kinetic equation (see TKB, Eqs. 3, 18) to a delta- 
function injection. In Appendix Bl we show that in this standard diffusion approximation 

< AE > is proportional to the velocity divergence, namely 

< A£ -1^<XZ£)b„ . Mb„ (11) 

3 K " 3r " ^ ' 

where k is the scattering coefficient (or the inverse of the scattering mean free path). But 
this method provides rather a qualitative estimate of the redshift effect in the outflow than 
a quantitative one. The value of the numerical factor 4/3 in equation (fTTl) is not precisely 
correct, which suggests that factor 4/3 should be replaced by the right one; however, the 
form of the exact equation actually depends on (3 = V/c and r in a more complex way. 
Formula ffTTj) should be then taken only as an estimate. 

A useful result of the < AE > estimate is obtained using the scattering geometry 
method (see Appendix B2) 



< AE >^ -/3//(2r)v/l-(//2r)2£;o. 



(12) 



Below we utilize the MC simulations to show that the numerical factor / in formula f|T2l) 
is about one. Comparison of Eqs. ffTTj) and f|T2l) leads us to conclude that < AE > is 
proportional to (3 but that is overestimated in equation Eq. flTT]) . 



The average energy of photons escaping after 
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1/6 scatterings is 



<E> 



sc 



(1+ < >) 




(13) 



From here using formula 0121] we obtain that 



< E {1 - /3//(2r)v/l-(//2r)2}V^Eo. 



(14) 



3.3.2. MC simulated spectral lines in outflow. Evidence of the strong redshift effect in the 



We simulate line profiles for different optical depths of the wind. In Figure 8 (left hand 
side panel) we present examples of the line profiles for simulations where we take into account 
only the pure electron scattering in the flow, photo-electric absorption being neglected in 
these simulations. In Figure 8 (right-hand side panel) we present the spectral profile of 
the line obtained as a result of the electron scattering in the flow along with photo-electric 
absorption. When the X-ray continuum spectrum of the central source illuminates the wind 
the seed i^^o— photons are formed in the innermost part of the wind. The photo- absorption 
profile as a function of energy, calculated by LT04, is related to the specific case of /3 = 0.1, 
a = 0.5, kTe = 50 keV for which the inferred LT04 model parameter TinT/L^Q = 6 x 10^^ 
cm. The model and X-ray continuum parameters are the following: r^^ is the outflow inner 
radius, L40 is the luminosity of the central source in units 10^" erg s~^, a is the energy 
spectral index and kT^. is the electron temperature of the Comptonization spectrum (see 
details of the photo-absorption calculations in LT04). We find that the height scale of the 
source exponential distribution, if, is only a small fraction of the inner radius or in other 
words H/rin ~ 0.04r is a small fraction of the optical thickness of the outflow. Thus the 
sources of is distributed in the innermost part of the flow. We simulate the line profiles 
for r = 2, 4, kT^ = 0.1 keV, (3 = 0.1, and for an initial seed photon energy Eph = 6.6 keV. 

One can clearly see the red-skewed part of the spectrum formed by photons undergoing 
multiple scatterings, along with the primary peak formed by photons escaping directly to 
the observer. The normalization of the direct non-scattering line component is suppressed 
for r = 2. The real normalizations are = 1.35 x 10^, 1.29 x 10^ for the diffusion and 
diffusion+absorption cases respectively. 



outflow 



In a situation, when the effect of the photon photoionization and absorption is taken into 
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account, the red wing of the hne (for energies less than 4 keV) and blue wing (for energies 
higher 8 keV) have to be suppressed due to high photo-electric absorption opacities. Despite 
this photo-absorption suppression, the red-skewed shape of the line is retained. Thus one 
can conclude that the red-skewed line profiles are intrinsic signatures of outward propagation 
of Ka— photons generated at the bottom of the outflow due to the outflow illumination by the 
X-ray hard radiation of the central source. 

In Figure 9 we show timing properties of the line propagation in the outflow for 5 energy 
bands for the pure scattering case. In the plot the time is given in the free path time units 
tfp = l/c. Photons escaping from the bottom atmosphere without scattering are at 6.6 keV. 
The model parameters used here are kTe = 0.1 keV, (3 = 0.1, r = 4 (upper panel), and 
r = 2 (lower panel). We show the arrival time of photons at the top of the cloud for five 
energy bands, in free path time units tfp. Soft lags, low energy emission at later times are 
clearly seen. 

It is also worth noting that TKB predict soft lags, namely that a photon of initial energy 
Eq loses its energy with time, i.e. E = exp(— 'Utkb)-E'o [compare with our formulas (IT^ 
]. The dimensionless time variable (we call it as mtkb) used there is a mathematical variable 
for a convolution method (see their Eq. 9). It can be interpreted as a product of the number 
of scatterings t/t/p multiplied by the mean photon inelasticity < AE > / E (see our Eq. [T2|) . 
i.e UTKB = (< > /E)(t/tfp). 

The photon distribution with respect to escape time integrated over energies can be 
fitted with an exponential law bexp{—bt/tfp) (shown in Fig. 2). The values of 6 = 0.057 
and b = 0.144 for r = 4 and r = 2 respectively (and consequently the average number of 
scatterings Nav = 1/fe) are very close to b (or Nav) obtained from the best-fit analytical 
spectrum for the same values of (3 and r. 

We also check the analytical estimate of the average redshift of the line energy using 
formula f|T^ and that obtained from our simulations with the values < E >sc= 4.56 keV 
and < E >sc= 5.27 keV for r = 4 and r = 2 respectively. On the other hand keeping in 
mind the values of b, related to a given r we can obtain the same values of < £" >sc as the 
simulated ones using formula f[T^ with / = 1.3, / = 1.7 for r = 2 and r = 4 respectively. 
These values of / > 1 are expected in the framework of our analytical estimate of < E >. 
(see details of the derivation in Eq. (IBllI) . It is worth noting that now we can calculate 
e =< AE > /Eq that follows from formula (fT2l) using / inferred from our simulations. They 
are e = 0.03, and 0.04 for r = 2, 4 respectively and for (3 = 0.1. 

In Figure 10 we present the emergent line profiles for different values of r and jS. Spectral 
lines as a function of r are for (3 = 0.1, and that as a function of (3 for r = 2. The source 
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photons are generated at the bottom of the cloud and the cloud temperature is kept fixed 
at kTf. = 0.1 kcV. The red wing of the line profile is insensitive to the temperature of the 
fiow if kTe < 1 keV. On the other hand, the blue wing becomes broader with kT^, but its 
relative width is still smaller than that of the red wing. In fact, the blue wing is affected 
for energies higher than K-edge (> 8 keV) if the photo- absorption is taken into account (see 
the left lower panels of Fig. 10). In general, the photo-absorption narrows the line profile 
(compare the lines profiles in the upper and lower panels of Fig. 10). 

The line profiles depend quite strongly on r and j3. The shape of the red wing below 
the broad peak (sec Fig. 10) follows a power law with index which is a strong function of 
r and j3. Our calculated indices are in agreement with the TKB results for (3 < 0.5 and 
Te// = 3/9t > 1. In TKB, the photon indices are equal to F = 2-|-A^, where are eigenvalues 
calculated using their equation (30). 

When (3^1 the spectrum of scattered photons becomes fiat even for r = 2. The 
relativistic spherical outflow completely washes out the bump induced by scatterings (see 
the right panels of Fig. 10). However, for relatively small /3, the power law tails still remain 
steep for higher r (see the left panels of Fig. 10). This is very similar to the results obtained 
by TKB. 

4. Compcirison of MC simulated line profiles with the data 

In Figures 11 and 12 (the left panel) we present the results of fitting our model to the 
XMM (EPIC) data of Wilms et al. (2001) obtained during the observation of the bright 
Seyfert 1 galaxy MCG-6-30-15. These observations were focused on the broad Fe Ka fine 
at ~ 6.4 keV. For MCG 6-30-15 we fit the continuum with a power law of index 1.8 and 
compare the residuals with our model. The best-fit parameters of the pure scattering model 
we found are Eph = 6.51 keV, kT^ = 0.1 keV, r = 1.2, f3 = 0.02. We generated the primary 
photons at the bottom of outfiow. If we take into account the LT04 model for the photo- 
absorption along with electron scattering, then the best-fit parameters become Ep^ = 6.51 
keV, kTe = 0.1 keV, r = 1.7, (3 = 0.02. 

The best-fit model spectrum is very similar to the steep spectrum presented in Fig. 10 
for f3 = 0.01 (black curve in the lower panel). The quality of the fits strongly depends on the 
outfiow optical depth r and bulk velocity /3. It is less sensitive to the outfiow temperature 
kTe unless kTgt 0.1 keV. It is important to note that using our diverging outfiow model we 
confirm Fabian's (1994) finding that velocity of the warm absorber outfiow is about 5000 km 
s-^ in MCG-6-30-15. 
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We fitted the red-skewed iron fine profiles for four galactic sources XTE J1550-564, GRO 
J1655-40, Cyg X-1 and GRS 1915+105 obtained by Miller et al. (2004b). AU of these data 
can be fitted by our model. In Figures 11 and 12 (the right hand panel) we show an example 
of these fits to the data. For GRO J1655-40 for the pure scattering model we find that optical 
depth r = 2.0 and an outfiow velocity is (3 = 0.1, much higher than that for MCG 6-30-15 
and Eph = 7.1 keV. If we take into account the photo-absorption along with scattering in 
our model, then the best-fit parameters are r = 3.3, /3 = 0.25 and Ep^ = 7.1 keV. It is worth 
noting that the fit in this case is not so good as that in the case of the pure scattering model. 
It can mean that the iron abundance in GRO J1655-40 is less than cosmic abundance. In 
fact, LT04 calculated the photoabsorption cross-sections for the cosmic abundances. 

As we have pointed out, the MG simulated line profiles are independent of the outfiow 
size and are determined by four parameters r, /3, kT^ and Eq. The size of the outfiow shell 
can be determined by timing characteristics such as variability timescales and time lags. Our 
model predicts soft time lags, i.e. the time lags become longer as energy decreases. Time lags 
are scaled with the light crossing time, tcross, from one to a few tcross- Unfortunately time lag 
information is not yet available for all sources analyzed. The only available information at 
present is the lack of the iron line variability with respect to the continuum for a number of 
extragactic sources [Markovitz, Edelson & Vaughan (2003)]. Thus one can suggest that these 
line features are formed in a more extended configuration than the source of the continuum 
(which is presumably a Gompton cloud close to the compact object). Here we demonstrate 
how our model can be used in the analysis of the line data. The detailed interpretation of 
the line data using our outfiow model will be presented in a separate publication. 

5. Outflow reprocessing of X-ray continuum spectra: MC simulated spectra 

We also study the modification of the central source X-ray spectrum in the warm outfiow 
as a result of photoionization and downscattering. In Figure 13 we present the simulations 
for different values of r = 2, 4 and /3 = 0.1, 0.3. The results of simulations are shown 
as E^F{E) diagrams, where F{E) is a photon spectrum. In these simulations, the incident 
(central source) spectrum is approximated by the sum of two components: a blackbody shape 
radiation presumably coming from the inner part of the accretion disk, and a power-law with 
an exponential cutoff coming from the Compton cloud (corona). The parameters used for 
this input spectrum are the color temperature of the blackbody component kThh = 1.2 keV, 
the power-law photon index F = 1.5 (energy index a = 0.5) , and the cutoff Ecutoff = 50 
keV. As shown in Gomptonization theory (see e.g. ST80, Titarchuk 1994) the exponential 
cutoff energy Ecutoff of the unsaturated Gomptonization spectrum is approximately 2kTe 
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where kTe is the plasma temperature of the Compton cloud. The temperature of spherically 
symmetric wind is kT^ = 0.1 keV. 

For the X-ray spectrum of the central source we assume that the photon numbers are 
the same for the blackbody and hard components. In fact, the ratio of the photon numbers 
depends on the illumination of the Compton cloud by the source of blackbody radiation. 

The photoelectric effect is taken into account in these simulations. When a photoelectric 
event occurred, we drew a random number to firstly determine if this event occurs in the K 
shell of iron. The probability of fluorescence is about 30 %. The fluorescence photon is then 
propagated into the cloud in the same way as the input central object spectrum. 

In the simulated spectra one can clearly see features of the strong fluorescent Kq, hne and 
the deep K-edge formed in the wind. The prominent bump around 25 keV is a combined 
effect of photoionization and downscattering in the wind. The wind modiflcation of the 

spectra is very sensitive to the optical depth of the wind r and (3. For example, for a given 
/3 = 0.1 (or the wind velocity 0.1c ) the depth of K-edge is broader and much deeper for 
r = 4 than for r = 2. As expected the width of the line increases with (3 (compare the 
spectrum for /3 = 0.1 with that for j3 = 0.3). It is important to emphasize that the shape of 
Ka line consists of narrow and broad components. These narrow and broad line features of 
the line are really observed in Cyg X-1 (see review of the line observations in Cyg X-1 
in the introduction section of ST06) . 

The shape of the spectrum for r = 4 and j3 = 0.3 is very similar to the observed 
spectrum of Cyg X-3 (TSh05 and Trudolyubov, 2006, private communication). We present 
the interpretation of Cyg X-3 data based on the results of our simulations elsewhere. 

6. Discussion and Conclusions 

It is important to point out that Thomson optical depth of outflow is an order of unity 
when the mass outflow rate Mout is an order of the Eddington mass accretion rate MEdd and 
higher (see e.g. formula 4 in King & Pounds 2003). 

Recent XMM-Newton observations of bright quasars (Pounds et al., 2003a,b; Reeves 
et al., 2003) give strong evidence for powerful outflows from the nucleus with mass rates 
Mout ~ Mq yr~^ ~ A^Edd and velocity V ~ 0.1c (i.e. {3 = 0.1) in the form of blueshifted 
X-ray absorption lines. These outflows closely resemble those recently inferred in a set of 
ultraluminous X-ray sources with extremely soft spectral components (Mukai et al. 2003; 
Fabbiano et al. 2003). It is also worth noting the recent observations of the strong outflows 
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from GRO J1655-40 (Miller et al. 2006) and the Nucleus of the Seyfert 1 Galaxy NGC 3783 
using an intensive HST/Chandra/FUSE monitoring (Crenshaw et al. 2004) and the recent 
Chandra observations of highly obscured AGN (Levenson et al. 2004) reported in the AAS 
HEAD meeting in 2004. In particular, Levenson et al. analyze new Chandra observations of 
three Compton thick Seyfert 2s with obscuring column densities that exceed 10^^ cm~^. They 
report that the Fe line is very prominent in these examples with equivalent width > 1 
keV. Thus one can conclude that powerful mass outflows from Eddington-limited accreting 
compact objects appear to be a very widespread phenomenon. King & Pounds (2003) and 
Begelman, King & Pringle (2006) argue that the powerful outflows may provide the high 
luminosity observed in quasars and ULXs, and imply that such objects have a major effect on 
their surroundings. They further suggest that the powerful outflow phenomena have barely 
been explored, and the field promises to be fruitful. 

It is out of the scope of this paper to study the consequences of the X-ray radiation 
for hydrodynamics of these outflows. Our goal is to elaborate the spectroscopic tools for 
probing the outflow physical conditions such as column density (the outflow optical depth) 
temperature and iron ionization charge state. Moreover the photoionization modeling of the 
outflow (see LT04) permits the determination of its distance from the central black hole. 

We have presented the Monte Carlo simulations of the radiative transfer of monochro- 
matic photons and continuum photons within a diverging outflow, a problem which is also 
of observational interest, in view of the observations of broad, redshifted Fe lines and down- 
scattering bumps (so called "reflection features") in cxtragalactic and galactic black hole 
candidate spectra. We also provide an analytic description of the simulated spectra which 
allows us to understand, in detail, the effect of the dowscattering modiflcation of the central 
source spectrum. In the form presented the derived formulae of the resulting spectra can be 
easily used as an efficient analytical tool for the spectral data analysis. We demonstrate that 
the strength of the downscattering bump is determined by the average number of scattering 
that photons undergo in the wind. From our simulations of the pure scattering wind model 
we find that the average number of scattering Nav = suffered by the photons in the wind 
of optical depth r^2 is ^r^ for f3tO-l (see §3 for details). On the other hand the average 
number of scatterings decreases when the velocity of the wind increases. A larger frac- 
tion of the photons are rather carried out by the flow than scattered there when velocities 
are higher than 0.1c. 

The outflow leads to a redshift of the photon energy producing a broad, red component 
(see Fig. 1 for illustration of this redshift effect) not unlike those of the Fe Ka lines observed. 
We demonstrate analytically and numerically (using Monte Carlo simulations) that the mean 
photon energy change per scattering in the outflow < AE > is the flrst order effect with 
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respect of F/c, namely < /S.E >oc —V/c. We further emphasize that the red-skewed hne 
features appear when the hne photons are generated very close to the bottom of the outflow 
due to the outflow illumination by the X-ray hard radiation of the central source. 

An important result of our simulations (see also TKB's results) is the natural suppression 
of a blue wing in the iron line feature without the need to invoke any specific geometric 
arrangement for the emission. The illumination of the outflow from inside by the X-ray 
radiation, is natural for a compact source geometry, where the X-ray radiation originates 
from the innermost part of the source (a Compton cloud along with a disk). The powerful 
wind presumably starts in outskirts of the soiut:c (sec Everett & Ballantyne 2004 for the 
wind model). In our model the red-wing photons have undergone numerous scatterings 
in an (effectively) optically thick medium, with an associated first order redshift in each 
scattering. It is not surprising that the blue wing, as a result of the second order {v/cY 
effect, is weak for kTet keV and (3 — V/ct^.l. 

This red-skewed line is a natural consequence of the first order Doppler effect in the 
presence of multiple scattering and photo- absorption events in the wind. It does not require 
any particularly fine tuned geometric arrangement. 

LT04 find the self-consistent temperature and ionization structure of the wind shell as a 
function of the parameter (radius/luminosity, rj^r/ L40 where r^n is the inner outfiow radius 
and L40 is X-ray luminosity in units of 10^*^ erg s~^) for a given Thomson optical depth 
of the shell. It is evident that the ionization parameter L^o/nor^ is constant through the 
shell (i.e. for r > ri„) if the velocity of the wind is constant through the fiow. Thus the 
LT04 solution allows one to determine the size of the shell base for a given luminosity of 
which the red-skewed line is observed. LT04 derive the range of parameter rin/L^Q where the 
outflow-produced stable solution exist. The allowed rin/ L4Q values are concentrated around 
IO12 _ iQis cm. 

Our model predicts the soft lags to he a function of E/Eq, j3 and r. The combination 
of the red-skewed features along with the soft time lags that strongly depends on /3, r and 
kTe are intrinsic signatures of any diverging outflow. 

The demonstrated application of our outflow model to data points out a potentially 
powerful spectral diagnostic for probes of the outflow-central object connection in Galactic 
and extragalactic BH sources. Particularly, analyzing MCG-6-30-15 data (sec §4) we infer 
that bulk outflow velocity is about 6000 km s~^. On the other hand the OVII edge in the 
warm absorber of MCG-6-30-15 indicates that the flow along line of sight is about 5000 ±1000 
km s~^ (Fabian et al. 1994). Thus using a new X-ray spectroscopic method [our diverging 
(outfiow) model] we confirm the early finding (that used UV spectroscopy) that velocity of 
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the warm absorber outflow of MCG-6-30-15 is about 5000 km s ^. 
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A. Downscattering spectrum. Derivation of Analytical Spectra 

A.l. Power-law as a spectrum of the primary photons 

Here we demonstrate how formula ([2]) can be analytically calculated using the steepest 
descent method. At first to illustrate the main idea of this derivation we consider the simplest 
case when e = and <f{z) = when 

iy(°)(z,n) = (1/z-n)"-^ (Al) 

and integral ([2]) is presented by equation (Q, where V{u) = bexp{—bu). If we change the 
variable v = {1/z — u)z under integral ([6]) then the integral can be rewritten as follows: 

^e{z) = exp{—b/z) / v''-^exp{bv/z)dv. (A2) 

z Jo 

Because the dimensionless frequency z = E/rrieC^ <^ 1 for the energy range of the interest 
[E <^ 511 keV) and btO.l there is a wide energy range where A = 6/2; is a big parameter. 
Then the integral in formula (IA2p can be calculated analytically using the steepest descent 
method. In order to do it one should expand the integrand function f"""^ in the Taylor's 
series 

v-^ = 1+^1^(1-.)+^^ - ^^i' - (1-.)^+...+ ~ ~ ~ (i-.r+... 

1! 2! n! 

(A3) 

near the point v = 1 where exp{bv/z) has a sharp maximum. It is worth noting that series 
(1A3I) converges for all f < 1 and thus the series presentation of the emergent spectrum (]A2[) 
is exact for any X = b/z: 

= ^eM-l-mUW + ^A(A) + <^'"^f'"> L.(A) + ... 

+ '^-°"'-°'-'"-°' l„(A) + ...]. (A4) 
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where ^ 

Ln{X)= / {1 ~ v)'' exp{Xv)dv for n = 0, 1, 2... (A5) 
Jo 

If we introduce the new variable u = {1 — v)/z in integral flASP then formula flA4p can be 
rewritten in the form 

:Fe{z) = Az--[Mo{z-') + ^-^^zM,{z-') + ^-^—^^^z'M,{z-') + ... 

+ ^^~"^^'~"^-^"~"^ ."M„(.-^) + ...], (A6) 
nl 

where Mn{z^^) is defined by equation ([8]). 



A. 2. Power-law with exponential cutoff as a spectrum of primary photons 

In the general case of the incident spectra and e > formula ([2]) for the emergent 
spectrum can be rewritten as 

b 

J^^{z,e) = - W^'\z,vUmax,s{z)]^M-^v)dv (A7) 

where A = bumax,s{z)^^ and W^'^\z,u) = e~^'^^^^[ijj^'^\z,u)](f[tp^^\z,u)] and u = vUmax,e{z)- 
Because the integrand has a sharp exponential maximum at f = (or -u = 0) one can expand 
the non-exponential function W^^\z,u) in series over v (or u). That leads to the following 
formula 

J^,{z,e) = ^[W^'\z,0)MoiUrnax,e) + Wi'^^'\z , 0) M,{Umax,e) / 1^ + 



+ Wi'^'^'\z,0)M2{Umax,e)/2\ + ... + Wi'^'^^\z , 0) Mn{Umax,e) / n\ + ...], (A8) 



X 



where Wu {z,u) is nth partial derivative of W^^\z,u) over u and the moments M„ 

(x = Umax,e) of V{u) aS folloWS 

M„(x) = -x"exp(-x) +nM„_i(x)/6 for n = I, 2 ... (A9) 

For example for n=0, 1, 2 we obtain 

Mo(x) = 1 - exp(-te), (AlO) 

Mi{x) = -xexp{-bx) + (l/6)Mo(x), (All) 

M2{x) = -x^exp{-bx) + (2/6)Mi(x). (A12) 
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Because 

n-1 

iyi^)(")(z,0) = Wi^\z,Q) n[lnW^i^^«(^,«)]U« = 0) (A13) 

and 

W^'\z,{)) = z^{z) (A14) 
the analytical formula ( lASp for the emergent spectrum can be presented as 



n=l 



n—1 



i=0 



)/n\ } . (A15) 



Thus the resulting spectrum J^^{z,e) is a product of the incident spectrum ip{z) and the 
downscattering transformation function (expression in the curl parenthesis in Eq. IA15I) that 
depends on the properties of the outflow, velocity, optical depth via the parameter b. 

In fact, the more convenient way to calculate the resulting spectrum J^,^{z,e) is to 
calculate the derivatives Wu ''^^\z,u) exphcitly using the factorized form of Wu^^^\z,u), 
namely 

Wjf^^'\z,u)=Wi'\z,uMz,u), (A16) 

where 

^z,u) = [\nzoip{zo)]u (A17) 

and 

(p{zo) = ZQ°'exp{-zo/z^). (A18) 
Using Eqs. (]A17I , lAlSP we obtain a formula for ^{z,u): 

/ M , 1 — a dzn 1 dzn 
^z,u) = [lnzocp{zo)]u = -4£/3 + ^ - (A19) 

Zq OU Zif clu 

Then one can calculate Wu'^^^\z,u) as follows 

ra-l 

W^i"^ = E^n-i^«""'"'^'^«^ n = 2,3,... (A20) 

1=0 

where = n{n — l)..(n — k — l)/k\ and Wu"^ = W. Here and below we omit superscript e 
in order to simplify the notation for Wu^^"\z, u) and its derivatives. We also introduce new 
notations W^^^ = W^'\z,u)u and = 



Because 



^ = (5/3);.o + ;^o' (A21) 
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we can write Eq. flA19|l as 



= -4e/3 + (1 -a)(e/3 + ^o) 



1 dzp 
du 



and find the derivatives of ^{z,u) as 



$W(z,«) = (l-a)|^-l5^ for n = l,2,3... 

The second and higher derivatives of Zq over u are 

d'^zp ^ (5/3)^ + 2zq — 
dv? du du ' 



du''' 



i=0 



for n = 3, 4... 



(A22) 



(A23) 



(A24) 



(A25) 



Now all derivatives wiT^ in Eq. (lA20l) can be calculated very easily because zq = z aX u = Q 
by definition. 
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B. The mean photon energy change in the bulk outflow. Analytical 

evaluations 

B.l. Diffusion method 

There is a standard procedure to estimate the energy gain and loss per scattering using 
diffusion approximation (Prasad et al. 1988 and Titarchuk, Mastihiadis & Kylafis 1997, 
Appendix D). Let us assume that in the initial moment t = the monochromatic line of 
energy Eq emits at any point of the scattering medium (outflow). The energy change of 
the line < AE > due to scattering can be found using the time-dependent diffusion kinetic 
equation: 

Oti 

— = LrU + L,yn (Bl) 
ou 

where n{r,v) is the photon occupation number, L^, Ly the space and energy operators 
respectively (see TKB, Eq. 3) and u = net is the time measured in units of mean time 
between scatterings (see Rybicki & Lightman 1979, ST80). Our goal is to find a small 
change of the radiation field background which at the initial moment is 

z'^n\u=o = z'^no = z5{z — zq). (B2) 

The energy change per scattering < Az >=< AE > /irteC^ can be found by multiplication 
of the kinetic equation (IBip by 2;^, substitution of equation (lB2p in the right hand side and 
integration over z. Thus we have 

A p 1 V(V/c) .dn , 

— =<Az>= / - ^ ' ' z'—dz, B3 

where 6z/6u is a photon energy change per scattering. 
Integration by parts gives 



< Az >= _1Z(XM / :,'Sn^dz. (B4) 



00 

3„ 







3 

Substituting z^uq = z6{z — Zq) (see Eq. ( IB2I) in this equation, we obtain 

< -I^E,. (B5) 

Because 

V(V/c) = ^ (B6) 

r 

and K,r = a^UeT t for the electron number density of the constant velocity wind we can 
rewrite equation flB5|) in the form 

< AE >^ -\^E,. (B7) 
3 r 
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B.2. Geometry of photon scattering in the outflow 

In this section we evaluate < AE > using equation ffTOl) and geometry of photon scat- 
tering in the outflow. Ultimately, we should calculate the mean of the scalar product of 
(Vi — V2) ■ n/c (n is a unit vector along the path of the photon scattered at the next point, 
see Fig. 1). 

The length of the vector 

|AV| = |Vi - V2I = 2Ksin0/2 ^ ?;sine for e<l (B8) 

where 6 is an angle between Vi and V2. Let / is the mean free path between two consecutive 
scatterings at ri and r2 (see Fig. 1) then 

If |ri — r2|/ri = Ar/r < 1 (ri r2 = r) and Ar < Z^/2r then we can estimate 1 — cos6 as 
follows 

l-cosO ^ {l/ry/2, (BIO) 

Thus we can flnd that 

sin e^l/r = l/{N,aTr) = rj^r.^r) ~ //r (Bll) 

when Ijr < \ and where the numerical factor The precise value of the factor / can be 
evaluated using our MC simulations. To derive formula fIBlip we use the electron number 
distribution of the constant velocity wind (see formula [1]). 

The cosine of the angle 7 between AV = (Vi — V2) and n is 

C0S7 = cos'?/'cos(7r — 6^/2) + sin-?/; sin(7r — ^^/2) cos((/) — (/Jq) (B12) 

where ■?/' is a zenith angle between vector n and the normal to vector Vi directed toward 
vector V2 and and are meridian (longitudinal) angles of n and ISV respectively. Then 
the mean value of cos 7 over half of sphere located to the left from the vector Vi is calculated 
as follows 

■j^ /'27r /'7r/2 -j^ 

< C0S7 >= — / d'P \ cos 7 sin ■i/'fi-?/' = — cos(6'/2). (B13) 
27r Jo h 2 

Thus using Eq. lITO!) and Eqs. (IBSHBlSl) one can relate the mean energy change per 
scattering, < AE >, for photons undergoing quite a few scatterings in the flow with /? as 
follows: 

< AE >^ -/3//(2r) Vl - if/2TyEo. (B14) 
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Fig. 1. — On the left side: Schematic diagram showing wind geometry. The relatively cold 
outflow (wind) originates at the inner radius. The optical depth of the wind in the Fe 
K continuum and the electron scattering optical depth of the wind are of order unity. A 
photon emitted near the inner boundary and subsequently scattered by an electron moving 
with velocity Vi, is received by an electron moving with velocity V2 as shown. The change 
in frequency is V2 = ^1 (1 + (Vi — V2) ■ n/c) where n is a unit vector along the path of the 
mv photon scattered at the next point. In a diverging flow (Vi — V2) • n/c < and photons 
are successively redshifted, until scattered to an observer at infinity. The color of photon 
path indicates the frequency shift in the rest frame of the receiver (electron or the Earth 
observer). On the right side: In a converging flow (Vi — V2) - n/c > and photons are 
blueshifted. 
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Fig. 2. — Timing properties of the outflow. The photon distribution over escape time 
V{t/tfp), where tfp = Z/c is a free path time. This simulation is made for the incident 
spectrum as a power-law with exponential cutoff for kTf, = 0.1 keV, r = 4 and (3 = 0.1. It is 
worth noting that the index b = 0.057 of the exponential distribution exp{—bt/tfp) (red line) 
is very close to the best-fit parameter values 0.06 ± 0.005 for the energy spectrum. We find 
that in all cases the best-fit parameter b obtained from the photon distribution over escape 
time is very close to b obtained from the best-fit analytical spectrum. 
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Fig. 3. — The outflow shell is radially illuminated be the central source radiation. Photon 
distribution over deflecting angles with respect to the radial direction of the incident photons. 
Upper panel: r = 1, middle panel: r = 2 and lower panel: r = 4 and (5 — 0.1. In fact, 
wc find in our simulations that the angle distribution of the scattering component weakly 
depends on r. But there is a big difference in the normalization of the direct (non-scattering) 
component Ajsi for which the ratio of Aj^{t = 4:)/An{t) = exp(— 4 + r). 
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Fig. 4. — Emerging spectrum of a compact object through an outflowing plasma. The 
incident spectrum inside the outflowing plasma is taken to be a power- law ^p{E) oc E~°' 
for which the energy index a = 0.5. Outflow parameters are fcTg = 0.1 keV, r = 4 and 
P = 0.05, 0.1, 0.3, 0.8. The photon emerging spectrum (cross points) is compared to the 
incident spectrum (dash line) and the analytical downscattering spectrum (solid line) (see 
Eq. ED. For (3 = 0.05, 0.1, 0.3, 0.8, the best-fit parameters are b = 0.07, 0.07, 0.11, 0.45 
respectively. For all cases we fix the e = 0. One can clearly see a bump in the MC simulated 
and analytical spectra, in the low energy band around 10 keV for all (3t0-3. Indeed, in 
the purely diffusive case that is shown in this plot, the photon number is conserved, so the 
downscattered high energy photons, which are removed from the high energy part of the 
incident spectrum, are detected at lower energy. 
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Fig. 5.— The same as in Fig. 4 for r = 2. For /? = 0.05, 0.1, 0.3, 0.8, the best-fit 
parameters are 6 = 0.13, 0.15, 0.21, 0.62 respectively. 
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Fig. 6. — Emerging spectrum of a compact object through an outflowing plasma. The 
incident spectrum inside the outflowing plasma is taken to be fcomp{E) oc E'" exp{—E / E^) 
(see Eq. |9]) where energy index a = 0.5 and high energy cutoff ai E = E^, = 130 keV. 
Outflow parameters are kT^, = 0.1 keV, r = 4. The photon emerging spectrum (cross 
points) is compared to the incident spectrum (dash line) and the analytical downscattering 
spectrum (solid line) (see Eq. IA8I) . For (3 = 0.1, 0.3, 0.5, 0.8, the best-fit parameters are 
b = 0.06, 0.10, 0.18, 0.38 respectively. The low energy downscattering bump effect is not 
seen in the emergent spectra. 
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Fig. 7. — The same as in Fig. 6 for r = 2. For l3 = 
are b = 0.14, 0.21, 0.30, 0.59 respectively. 



0.1, 0.3, 0.5, 0.8, the best-fit parameters 
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Fig. 8. — Line profiles as a function of r (for jS = 0.1). Primary photons are generated 
at the bottom of outflow. On the left-hand side : the line profiles for the pure scattering 
case and on the right-hand side : that for scattering plus absorption case. In each case, 
the fixed parameters are Eph = 6.6 keV, kT^ = 0.1 keV. The normalization of the direct 
non-scattering line component is suppressed for r = 2. The real normalizations are An = 
1.35 X 10^, 1.29 X 10^ for the diffusion and diffusion+absorption cases respectively. 
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Fig. 9. — Timing properties of the line propagation in the outflow. A flash of monochro- 
matic photons at 6.6 keV is simulated, and propagated. Photons escaping from the bottom 
atmosphere without scattering are at 6.6 keV (see Fig. 1). The model parameters used here 
are kT^ = 0.1 keV, /3 = 0.1 and r = 4 (upper panel), and r = 2 (lower panel). We show 
the arrival time of photons at the top of the cloud for five energy bands, in free path time 
units tfp = l/c. Soft lags, low energy emission at later times are clearly seen. The photon 
distribution over escape time integrated over energies can be fitted with an exponential law 
bexp{—bt/tfp) (shown in Fig. 2). The values of b = 0.057 and b = 0.144 for r = 4 and r = 2 
respectively (and consequently average number of scatterings Nav = 1/&) are very close to b 
(or Nav) obtained from the best-fit analytical spectrum for the same values of (3 and r (see 
text). 
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Fig. 10. — Line profiles as a function of r and (3. Primary photons are generated in the 
bottom of outflow. Upper panel: the line profiles for the pure scattering case. Lower panel: 
that for scattering plus absorption case. In each case, the fixed parameters are Eph = 6.6 
keV and fcTg = 0.1 keV. Spectral lines as a function of r are for f3 = 0.1 and that as a 
function of /? for r = 2. Note, the case with a notation r = 0.0 [see upper and lower panels 
(a)], corresponds to r = 0.01. 



-37- 



0.06 



> 

OJ 



o 

-t— ' 

o 

CL 



0.04 



0.02 



0.00 - 



-0.02 





-0.2 



keV 



Fig. 11.— Fits of our model to the XMM observations of MCG 6-30-15 (Wilms et al. 2003) 
and to the ASCA observations of GRO J1655-40 (Miller et al. 2004b). For MCG 6-30-15, 
we have fitted the continuum with a power law of index 1.8 and compare the residuals with 
pure scattering model. The pure scattering model. For this model the best-fit parameters 
which we found are Eph = 6.51 keV, kT^ = 0.1 keV, r = 1.2, f] = 0.02. Primary photons are 
generated at the bottom of outflow (illumination of outflow from inside). For GRO J1655-40 
we found the following best-fit parameters: Eph = 7.1 keV, kT^ = 0.1 keV, r = 2.0, /5 = 0.1. 
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Fig. 12.— Fits of our model to the XMM observations of MCG 6-30-15 and to the ASCA 
observations of GRO J1655-40 using the model for which the photo-absorption along with 
scattering are included. For MCG 6-30-15 the best-fit model parameters which we found 
are Eph = 6.51 keV, kT^ = 0.1 keV, r = 1.7, /? = 0.02. For GRO J1655-40 we found the 
following best-fit parameters: Eph = 7.1 keV, kT^ = 0.1 keV, r = 3.3, j3 = 0.25. 
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Fig. 13. — The central source spectrum reprocessed through the wind shelL Upper panel: 
(3 = 0.1 and tq = 2 and tq = 4. Lower panel: that is for (3 = 0.3. In the X-ray spectrum of 
the central source the photon numbers are the same in the blackbody and hard components. 
In fact, the ratio of the photon numbers is arbitrary and it depends on the illumination of 
the Compton cloud by the source of blackbody radiation. In this simulations we assume the 
blackbody color temperature is 1.2 keV. The resulting spectrum is shown by black histogram, 
whereas orange and pink curves present the hard and blackbody components of the central 
source reprocessed in the wind respectively. A blue curve is the Kq line formed in the 
wind. One can clearly see the strong fluorescent Kq, line, K-edge formed in the wind and 
a prominent bump around 25 keV. The line consists of narrow and broad components (see 
narrow blue peak at 6.4 keV in all spectra). 



